Optimal spin-entangled electron-hole pair pump 
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A nonperturbative theory is presented for the creation by an oscillating potential of spin-entangled 
electron- hole pairs in the Fermi sea. In the weak potential limit, considered earlier by Samuelsson 
and Biittiker, the entanglement production is much less than one bit per cycle. We demonstrate 
that a strong potential oscillation can produce an average of one Bell pair per two cycles, making it 
an efficient source of entangled flying qubits. 
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The quantum electron pump is a device that transfers 
electrons phase coherently between two reservoirs at the 
same voltage, by means of a slowly oscillating voltage 
on a gate electrode pj. Special pump cycles exist that 
transfer the charge in a quantized fashion, one e per cy- 
cle 

HHHH Building on earlier proposals to stochasti- 
cally produce entangled electron-hole pairs in a Fermi sea 
out of equilibrium 0, , Samuelsson and Biittiker have 
proposed 8] that a quantum pump could be used to cre- 
ate entangled Bell pairs in a controlled manner, clocked 
by the gate voltage. Such a device could be a building 
block of quantum computing designs using ballistic fly- 
in g qu bits in nanowires or in quantum Hall edge channels 

To find out how close one get to this ideal, one needs to 
go beyond the perturbation theory of Ref. |8| — in which 
the number of Bell pairs per cycle is <C 1. A nonpertur- 
bative theory of the quantum entanglement pump is pre- 
sented here. We show that the entanglement production 
is closely related to the charge noise, to the extent that a 
noiseless pump produces no entanglement. By maximiz- 
ing the charge noise with spin-independent scattering we 
calculate that a pump can produce, on average, 1 Bell 
pair every 2 cycles. A deterministic spin-entangler [l^ . 
being the analogue of a quantized charge pump, would 
have an entanglement production of 1 Bell pair per cy- 
cle, so the optimal entanglement pump has one half the 
efficiency of a deterministic entangler. 

We consider a two-channel phase coherent conductor, 
see Fig.^ connecting a left and a right electron reservoir 
in thermal equilibrium (same temperature T and Fermi 
energy Ep in each reservoir). The two channels may refer 
to an orbital or to a spin degree of freedom. (To be defi- 
nite, we will usually speak of a spin degree of freedom.) 
A periodically varying time-dependent electrical poten- 
tial V(r, t) (with period 2-k/uj) excites electron-hole pairs 
in the Fermi sea of the conductor. The quantum me- 
chanical state of an electron-hole pair, at energies E, E' 
differing by a multiple of Tito, may be entangled in the 
channel indices. The entanglement can be a resource for 
quantum computing if the electron and the hole excita- 
tion are scattered to opposite ends of the conductor, so 
that they become two separate qubits. We wish to relate 
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FIG. 1: Production of entangled electron-hole pairs in a nar- 
row ballistic conductor by a quantum electron pump. The 
left and right ends of the conductor are at the same poten- 
tial, while the potential on the gate electrodes at the center 
is periodically modulated. Such a device produces spatially 
separated electron-hole pairs (black and white circles) , differ- 
ing in energy by a multiple of the pump frequency to. For 
spin-independent scattering the electron (e) and hole (h) pro- 
duced during a given cycle have the same spin 1,1, so that 
their wave function is that of a Bell pair, oc | tetft) + I ieifc)- 
The optimal quantum entanglement pump produces, on av- 
erage, 1 Bell pair every 2 cycles. 



this entanglement production to the scattering matrix of 
the conductor. 

The four characteristic energy scales of this problem 
are the thermal energy ksT, the pump energy fiuj, the 
Thouless energy Ti/rr> (set by the inverse of the mean 
dwell time Tr> of an electron in the conductor), and fi- 
nally the Fermi energy Ep. In nanostructures at low tem- 
peratures, the characteristic relative magnitude of these 
energy scales is ksT -C hid <C Ti/td -C Ep. This is the 
adiabatic, low temperature regime in which we will work. 

We seek a relation between the entanglement produc- 
tion and the scattering matrix S of the pump, which is 
the unitary operator relating incoming to outgoing states: 
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Here a n (E) is the fermion annihilation operator for an 
incoming channel n at energy E, and b n (E) is the an- 
nihilation operator for an outgoing channel. There are 
four channels in total (n — 1,2,3,4), two in the left lead 
and two in the right lead. The Wigner transform of the 
scattering matrix [l^, defined by 
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depends on E on the scale H/td- In the adiabatic regime 
lutjj <C 1 one may therefore neglect the ^-dependence on 
the scale of the pump energy, approximating Sw(E, t) ~ 
Sw{EF,t) = S(t). The 4x4 unitary matrix S(t) can be 
obtained by solving the static scattering problem for the 
frozen potential V(r, t). 

Within a single period 2tt/uj the excitation energies 
can only be resolved on the scale of hto, so we discretize 
E p = phuj with integer p. A pair of energies E p and E q 
is coupled by the Floquet matrix T(j> — q) , which is the 
Fourier transform of the Wigner transformed scattering 
matrix [l4| . 

Hp -i) = ^ J 2 *" dtS{t)e^-^K (3) 

The unitarity relation for the Floquet matrix reads 

E T nn> (P ~ P')Kin> (<? ~ P') = <>n m S pq . (4) 
n' ,p' 

We assume zero temperature, so the incoming state 
|^o) is the unperturbed Fermi sea, consisting of all levels 
(left and right) doubly occupied below Ep and empty 
above Ep: 



\*o)=l[f(E p )l[at(E p )\0). 



(5) 



The state |0) is the vacuum and f(E) — 6(Ep — E) is the 
zero-temperature Fermi function. The outgoing state | '&) 
is obtained from |$ ) by substituting Eq. JIJ and taking 
the adiabatic approximation for the scattering matrix, 

i $ > = II^)II(EE b l ( e p' )H>n( P ' -p))\o). 

p n ^ p' n' 

(6) 

We denote by w pq the probability that the pump ex- 
cites within one cycle a single electron-hole pair, consist- 
ing of an electron at the left at an energy E p above the 
Fermi level and a hole at the right at an energy E q below 
the Fermi level. The entanglement entropy (or entangle- 
ment of formation) of the spins is denoted by £ pq (mea- 
sured in bits per cycle). Similarly, w qp and £ qp refer to 



a hole at the left and an electron at the right. The aver- 
age production, per cycle, of spin-entangled electron-hole 
pairs is 
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A maximally entangled Bell pair has £ pq — 1, so it con- 
tributes w pq bits to the entanglement production. 

The weight factor w pq and entanglement entropy £ pq 
of an electron-hole pair can both be calculated by pro- 
jecting l^) onto a state V p g\^) with all levels (left and 
right) empty above Ep and doubly occupied below Ep 
— except for a singly occupied level E p > Ep at the left 
and E q < Ep at the right. The (unnormalized) projected 
electron-hole state has the form 

n»= E a ^ b L( E P ) b L'(Eq)\o). (8) 

The four channels have been labeled L j , L j , B j , R ! J. , 
where L, R refers to the left and right lead and the arrows 
t, | indicate the spin. The 2x2 matrix a determines the 
weight factor as well as the entanglement entropy, 
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(9) 

T ; ( 10 ) 
Ir act* 

The number C € [0,1] is the concurrence of the 
electron- hole pair. 

In order to calculate the matrix a it is more convenient 
to perform the algebraic manipulations on the pair cor- 
relator K in the outgoing state |$), rather than on the 
state itself. The pair correlator fully characterizes the 
outgoing state 10 because it is Gaussian, meaning that 
higher order correlators in normal order (all 6^'s to the 
left of the b's) are constructed from the pair correlator ac- 
cording to the rule of Gaussian averages. The correlator 
is given in terms of the Floquet matrix by 

K nm (p,q) = (bl(E q )b n (E p )) 

= J2 :F n,Ap-p')f(E p ')KirA<l-p')- (11) 

n' ,p' 

The matrix K is Hermitian and idempotent in the joint 
set of energy and channel indices: K = = K 2 . This 
signifies that the state it represents is a pure (rather than 
a mixed) state [l6j |. 

Projection of \^>) onto a set of filled or empty levels 
preserves the Gaussian property. The correlator K of 
the projected state VpqY^) is derived from K by the pro- 
cedure known in matrix algebra as Gaussian elimination 
[l7j . By interchanging rows and columns of the matrix 
K we move the indices pL^,pLl, qR\, qR[ to the upper 
left hand corner, to obtain the block form 



K 



At 



fK LL {p,p) K LR (p,qy 
\K RL {q,p) K RR (q,q) 



(12) 
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The 4x4 block K^h- contains the direct coupling of the 
spin degenerate levels p at the left and q at the right. The 
correlator K contains in addition the indirect coupling 
via the filled or empty states in the block Y, 

K = K dil + X(l-Y-A)- 1 X\ (13a) 
«>=(*|7' pg |*) = |Det(l-y-A)|. (13b) 

The diagonal matrix A has a 1 on the diagonal if the 
state is filled (below Ep) and a if it is empty (above 
Ef). A derivation of Eq. (|13|) is given in Add. 1X1 

One readily verifies that K 2 = K, so the projec- 
tion preserves the purity of the state, as it should. 
Since V pq \^>) contains a total of two electrons in four 
states, the correlator K has two eigenvalues equal to 
1 and two eigenvalues equal to 0. We write K — 
f/diag (0, 0, 1, 1)U>, with U the unitary matrix of eigen- 
vectors. The projected state corresponding to the corre- 
lator K has the form 

V pq \^>) = w^ 2 (b^U) m (^U) Rl \0), (14) 
b= (b L1 (E p ),b Ll (E p ),b R1 (E q ),b Rl (E q )). (15) 

The matrix U plays the role of an effective scattering ma- 
trix for the two electrons in the two states left and right, 
including in addition to the direct scattering (described 
by the original scattering matrix S) also the indirect tran- 
sitions via the other states. 

To obtain the required projection V pq \^) we still need 
to project V pq \^) onto a state with a single electron left 
and a single hole right, excluding the double occupation. 
(We could not do the projection in a single step because 
the final state (jHJ is not Gaussian, so it can not be repre- 
sented by a pair correlator.) By comparing Eqs. JSJ and 
(|14|l we can relate the coefficient matrices U and a before 
and after projection, 

a = i^U LR a v Ul R , U=(°£ (16) 

(The matrix a y is a Pauli matrix.) Substitution into 
Eqs. 10 and 11(J|I then gives the contribution from this 
electron-hole pair to the entanglement production. 

A major simplification occurs in the case of spin- 
independent scattering. Then Ul R and U RR are pro- 
portional to the 2x2 unit matrix 1, so a oc a y and the 
electron- hole pair is maximally entangled (£ pq =1). In 
view of Eq. Q the average entanglement production per 
cycle, 

e = E E «+<)> (it) 

is the probability that the pump produces a single spa- 
tially separated electron-hole pair in a given cycle. 

The probability (|17fl can be rewritten as £ = PqP\ + 
PqPi, where P° is the probability that v spatially sepa- 
rated electron-hole pairs of spin a are produced in a given 



cycle. From < P{ = P{ < 1 - Pj = 1 - P£ < 1 we 
deduce that 

£<2P CT (1-P CT )<±. (18) 

This maximal entanglement £ max = 5 of one half bit per 
cycle is reached for P T = P X = P/ = P/ = |. Eq. 
(|18f) is derived for spin-independent scattering. It seems 
unlikely that spin-dependent scattering (which reduces 
the entanglement per electron-hole pair) could violate the 
bound £ < i, but we have not been able to exclude this 
possibility on mathematical grounds. 

To demonstrate that the optimal value £ max = \ can 
be reached, we consider the pump cycle 

Sir) = e -i* V ) , (19) 

which has been used as a model for a quantized charge 
pump 0, E| • (A more general class of pump cycles 
with the same entanglement production is analyzed in 
App.IbI) Choosing the Fermi level such that Eq < Ep < 
Ei , Eq. ifTTjl evaluates to 

K{p,q) - (I f)5 pq f(E p ) 

f-rrUpoSgo tr n 5 p0 S ql \ , . 
+ \r'tU pl 5 q0 r'r'tvW' 

The only pair of coupled levels is Eq at the left and E\ 
at the right, so the entanglement production consists of a 
single term £ — w^fE^f. The matrix X in the decompo- 
sition (|12(l vanishes, and Y = A, hence Eq. (|13fl simplifies 
to 

K = K A „=(^-f w = l. (21) 

Eq. (|16f) gives a = ita y r' T , which finally leads to the 
entanglement production 

£=H(x 1 ,x 2 ), x 1 =T 1 (l-T 2 ), x 2 = T 2 (l-T 1 ), (22) 

in terms of the function 

H(x,y) = (x + y)\og 2 {x + y) - x\og 2 x - ylog 2 y (23) 

of the two transmission eigenvalues T± , T 2 (eigenvalues of 
tv ', equal to the eigenvalues of t't'> because of unitarity 
of S). 

The optimal entanglement production £ mayi = \ is 
reached in Eq. (|22|l for T\ = T 2 — \ (corresponding 
to spin- independent scattering, as expected). This is 
also the choice of parameters at which the charge noise 
oc Ti(l - T x ) + T 2 (l - T 2 ) is maximized 18]. Although 
entanglement entropy and charge noise are different phys- 
ical quantities, with a different dependence on the trans- 
mission eigenvalues, quite generally one can state that 
there can be no entanglement production without charge 
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noise. Indeed, a deterministic spin-independent charge 
pump has Pq = hence £ = 0, in view of Eq. i|18fl . 

A one-to-one relationship between entanglement pro- 
duction and charge noise is possible in the weak pumping 
limit of Ref. To demonstrate this, we quantify the 
pumping strength by a dimcnsionless parameter e <C 1, 
and calculate both quantities to leading order in e. The 
Floquet matrix to first order has the general form 

To if p = q, 

Hp -q) = { ^QH if v = q + 1, (24) 

ieQ^J-o if p = q — 1. 

Unitarity of J-q ensures unitarity of T up to terms of 
order e 2 . The corresponding correlator Ulljl is 

K(p,q) = f(E p )S pq + ie[f(E q ) - f{E p )]5 p , q+l Q 

- ie[f(E p ) - f(E q )]S q , p+1 Ql (25) 

Following the same steps as before, we arrive at the 
entanglement production 



£ = e 2 H(y l7 y 2 )+e 2 H(y[,y' 2 ) 



(26) 



in terms of the function H defined in Eq. J2HJ). The 
numbers y n and y' n are the eigenvalues of the matrices 
tt^ and t't'\ respectively, constructed from sub-blocks 
of the matrix 



P T 



(27) 



In the case of spin-independent scattering yi = y 2 



y', and Eq. ffify simplifies to 
£ = 2e 2 (y + y' 



(28) 



To compare this result with the charge noise V, we use 
the formula H S E3 



OO 

P = -^pTrG(p)G(-p), 



(29) 



G[p) 



E -*>)(! °t)Hq)- (30) 



q— — oo 



In the weak pumping limit this reduces to 



limit to the frequency. What motivates further efforts 
on the side of the detection is the relative simplicity on 
the side of the production: The quantum entanglement 
pump requires no advanced lithography or control over 
electron-electron interactions to produce as much as one 
Bell pair per two cycles. 

This work was supported by the Dutch Science Foun- 
dation NWO/FOM and by the U.S. Army Research Of- 
fice (Grant No. DAAD 19-02-0086). 



APPENDIX A: GAUSSIAN ELIMINATION OF 
FERMION DEGREES OF FREEDOM 

The Gaussian elimination of boson degrees of freedom 
was done in Ref. 0] . Here we perform the analogous cal- 
culation for fermions, leading to Eq. I|13f) . We represent 
the density matrix p by the normally ordered character- 
istic function 



x(0 



Irpe q e ? 



(Al) 



The entire set of fermion annihilation operators is con- 
tained in the vector b. The characteristic function de- 
pends on the vector of Grassmann variables For a 
Gaussian state, with pair correlator K, one has 



X(0 = exp(£* 



(A2) 



We first eliminate a single degree freedom, labeled by a 
subscript 0, by projecting onto a filled state. The projec- 
tion operator is &J& , so the characteristic function after 
projection is 



X0,filled(£) 



= Trblb Q pe^ b1 e^ b 

V d d ft\ 

to, iS^o d£ 



dHo / d£Sx(€)- 



(A3) 



The final equality holds because differentiation and inte- 
gration is the same for Grassmann variables. Similarly, 
we can project onto an empty state, with projection op- 
erator 1 — b b and characteristic function 



V = e 2 (y 1 +y 2 + y' 1 +y' 2 ) 



(31) 



which equals the entanglement production (|28fl in the 
spin-independent case. 

The close relation between entanglement production 
and charge noise in the weak pumping regime is consis- 
tent with the finding of Ref. [8| that low-frequency noise 
measurements can be used to detect the entanglement in 
this regime. To access the nonperturbative regime inves- 
tigated in this paper requires time-resolved detection, on 
the time scale of 1/w. The requirement that the thermal 
energy ksT remains less than Tiuj poses a practical lower 



X0,ompty(€) = Tr (1 - b\b Q )p e« bt e«* b 



■fdto /d8e-«*o x (£). (A4) 



We now divide the degrees of freedom into two sets, 
those that are to be eliminated by projecting onto filled 
or empty states (labeled P), and those that are to be re- 
tained (labeled R). The diagonal matrix A, in the space 
of projected states, has a on the diagonal if the state 
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is empty and a 1 if the state is filled. The characteristic 
function after projection is 



x e 



-^•(l-A).^ 



x(0- 



(A5) 



The Gaussian characteristic function l|A2|) has the 
block structure 



X(£n,£p) = exp 



Vp 



K dil X 
Xt Y 



Sit 
Zp 



(A6) 



Substitution into Eq. I)A5|I gives, upon integration over 



J 



the Grassman variables, 



= II (-l) 1 "^" Det (V — 1 + A) 



riGP 



x exp 



e R ■ (K dil -X-(Y-1 + A)- 1 • Xt) • & 



(A7) 



This is the characteristic function of a Gaussian state 
with weight and pair correlator given by Eq. {15} ■ Wc 
note that this derivation holds regardless of whether 
p represents a pure or a mixed state (so regardless of 
whether K 2 = K or not). 



APPENDIX B: CLASS OF OPTIMAL ENTANGLEMENT PUMP CYCLES 



We calculate the entanglement production for pump cycles of the form 



s(r)=r:r ^ )= j_^_ ( bd 



t' e-^r'J ' 1 - z*e 41 

introduced by Ivanov, Lee, and Levitov 19] in the context of charge pumps. The complex parameter z satisfies 
\z\ 2 < 1. We will show that these pump cycles all have the same entanglement production as the special case z = 
considered in the main text [Eqs. (|19|) - H23|l ]. In particular, they are optimal for any z if Ti = T% = h. 

We label the energy levels such that E p < Ep if p < and E p > Ep if p > 1. The pair correlator is given by 

K(v <,)-(* °} 6 f + (1 \z\ 2 ) ( -^f P f q z- q (z*r p try p (l - f q ){z*Y-^ \ 

We have abbreviated f p = f(E p ). This pump transfers electrons exclusively from left to right, so only half of the 
terms in the entanglement production J7J) need to be included: 

oo 

p— — OO <J— 1 

The index p < refers to a hole excitation at the left and the index q > 1 to an electron excitation at the right. 

The Gaussian elimination is more complicated than for the case z — considered in the main text, because 
the coupling between two given levels p and q now contains also an indirect contribution (via an infinite ladder of hole 
excitations at the left and electron excitations at the right). Still, after some algebra the final result for the correlator 
K of the projected state V pq \^) turns out to be rather simple. We introduce the polar decomposition 



r t \ _ (u 0\ (VT^T -VT \ fu' 
f r'J - [O v) I VT VT^TJ ( v' 



(B4) 



where u,v,u',v' are 2x2 unitary matrices and T = diag(Ti,T2) is a diagonal matrix of transmission eigenvalues. 
Defining also c pq = (1 — \z\ 2 )z q ~ p ~ 1 , the result for K can be written in the form 



<">)\ T+(\-T)\c„\") - T) (l-T)M 2 J\0 »' 

The weight factor w equals 

w= 11 [T n + (1 - T n )\c pq \ 2 ] . (B6) 



71=1,2 
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The unitary matrix U that diagonalizes K — [/diag (0, 0, 1, 1)W is given by 



u 0\ {T+(l-T)\c pq \ 2 \ 1/2 fc*VT=T VT 



" = r :) [' " , r+(1 : TM ) ^ ^fcr, . (B7) 



hence Eq. iJTSj l gives 



a = ^ /2 ° pqU ( t + (i T T)| Cp9 |0 7 ff * (r + (i-T)| Cp9 p) 7 vT - (B8) 



Substitution of Eq. (|B8|I into Eqs. © and {TU1) gives 
the weighted entanglement production of the electron- 
hole pair, 

<£p 9 e = \c Pq \ 2 H (Ti(l - T a ), T a (l - Ti)), (B9) 

in terms of the function H defined in Eq. i|23|) . Finally, 
we carry out the sum over the level indices p and q in Eq. 



(|B3|) to obtain the total entanglement production 

£ = H(T X {1 - T 2 ),T 2 (1 - Ti)), (BIO) 

which is the same as Eq. (|22H . The z-dependence drops 
out because X]p=-oo S^Li \ c pq\ 2 = 1- The maximum 
£ = h is reached for Ti = T 2 = 5, independent of z. 
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